Mon. Not. R. Astron. Soc. 000,[T]tT2](201 1) Printed 4 July 201 1 (MN KT^ style file v2.2) 

Varying disc-magnetosphere coupling as the origin of pulse profile 
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ABSTRACT 

Accreting millisecond pulsars show significant variability of their pulse profiles, especially at 
low accretion rates. On the other hand, their X-ray spectra are remarkably similar with not 
much variability over the course of the outbursts. For the first time, we have discovered that 
during the 2008 outburst of SAX J 1 808.4-3658 a major pulse profile change was accompanied 
by a dramatic variation of the disc luminosity at almost constant total luminosity. We argue 
that this phenomenon is related to a change in the coupling between the neutron star magnetic 
field and the accretion disc. The varying size of the pulsar magnetosphere can influence the 
accretion curtain geometry and affect the shape and the size of the hotspots. Using this physi- 
cal picture, we develop a self-consistent model that successfully describes simultaneously the 
pulse profile variation as well as the spectral transition. Our findings are particularly impor- 
tant for testing the theories of accretion onto magnetized neutron stars, better understanding 
of the accretion geometry as well as the physics of disc-magnetosphere coupling. The identi- 
fication that varying hotspot size can lead to pulse profile changes has profound implications 
for determination of the neutron star masses and radii. 

Key words: accretion, accretion discs - pulsars: individual: SAX J1808.4-3658 - stars: neu- 
tron - X-rays: binaries 



1 INTRODUCTION 

SAX J1808.4-3658 was first detected with the Wide Field Cam- 
eras on board the BeppoSAX satellite in 1996 din 't Zand et al.l 
1 1998b . In 1998, another outburst of SAX J1808.4-3658 was ob- 
served with Rossi X-ray Timing Explorer (RXTE) and the discov- 
ery of ~ 401 Hz pulsations led to its identific ation as the first ac- 
creting millisecond pulsar (A MP; Iwiinands & van der Klijl 19981 ; 
IChakrabartv & Morganlll993) . Since 1998, the source has been in 
outburst roughly every 2.5 years (2000, 2002, 2005 and 2008). 
During the outbursts of SAX J1808.4-3658 (and other AMPs), 
the magnetic field of the neutron star (NS) channels the accreted 
matter on to the stellar magnetic poles. The emitted radiation 
from these hotspots is then modulated at the stellar spin period, 
that results in coherent pul sations. In a typical outburst of SAX 
J1808.4-3658 (see fig.2 in iHartman et al.ll2008 ') the source rises 
from quiescence in roughly five days, reaching a peak flux of about 

2 X 10"^ erg cm"^ s"^ in the 2-25 keV range of the RXTE Pro- 
portional Counter Array (PCA). After the peak, the source flux 
drops slowly ("slow decay" stage) in the course of the next ten 
days, although the 2008 outburst had a double peaked light curve 
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JHartman et al J 120091) . The slow decay stage is followed by the 
"rapid drop" stage, where the flux drops by an order of magni- 
tude in a couple of days after which the source exhibits several 
re-brightening episodes every five days or so ("flaring tail" stage, 
see also Patruno et al. 2009b, Ibragimov & Poutanen 2009). 

The system parameters of SAX 11808.4-3658 are the 
best known among AMPs. The orbital period i s ~ 2 hours 
l! Chakrabartv & Morgan|[T99i ; iBurderi et al1l2009l ; iHartman et all 
2009) an d the companion star has a very low mass of Mc ~ 
0.05Mq (Bildsten & Chakrabartv 20oT). The analysis of type I X- 
ray bursts by Galloway & Cumming ( 2006) gave a distance esti- 
mate of 3.5 ± 0.1 kpc. The inclination has been constrain ed to 
i =36°-67° using optical observations JDelove et al.l 120081) and 
the X-ray analysis jlbragimov & Poutanenll2009n of the 2002 out- 
burst led to a similar constraint of i « 50° -70°. By study i ng the 
long term timing evolution of the system, IHartman et al.l ( 1200 8h 
constrained the neutron star surface magnetic dipole field to a 
range of Bs = (0.4-1.5) x 10** G. Another estimate of Bs = 
(0.8 ± 0.5) X lO'^fcA G from pulse profile variability, was de- 
rived by Ibragimov & Pouta nen (2009), where the factor ^a ~ 0.5 
dLong et al.ll2005l) . Also. ISurderi et all (l200d) gave an e stimate of 
Bs ~ (3.5 ± 0.5) X 10* G and IPatruno et al.l ( l2009cl) obtained 
Bs = (2-3) X 10** G. 
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The energy spe ctra of SAX J1808.4-3658 has been stud- 
ied extensively (e.g. Gilfanov et al. 19981: Gierlinski et al.l 120021 : 
IPoutanen & Gierliriski 2003; Ibraeimov & Poutanen 2009). It was 
pointed out already by Gilfanov et al. ( 1998) that the spectral shape 
remains remarkably similar when flux changes by more than an or- 
der of magnitude throughout the outburst. The energy spectrum is 
roughly flat (photon index F ~ 2) with a cutoff/roll-over at en- 
ergies above ~ 50 keV. Such spectra - also seen in other AMPs 
JPoutanea uOOq) - are most likely produced by thermal Comp- 
tonization of soft seed photons originating from the heated stel- 
lar surface by hot electron s in the accretio n shock jG ierlihski e t al.l 
I2OO2I : IPoutanen & GierUh ski 2003; Ibrag imov & Po utanen 200^ 
Below ~ 5 keV an additional T ~ 0.5 keV thermal component is 
visible that is associated with the heated surface of the NS. XMM- 
Newton observations have revealed another non-pulsatine coole r 
thermal component due to the accretion disc fPatrun o et alj2009ah . 
The X-ray emission from the NS surface irradiates the accretion 
disc, that caus es spectral hardening above 10 keV due to Comp- 
ton reflection jlbragimov & PoutanenI 12009,) . This irradiation also 
produces an iron line at 6.4 keV and the modelling o f the line pro- 
files fr om XMM-Newton and Suzaku observations bv lPapitto et al.l 
( l2009f) and Cackett et al. (2009) has given tight constraints on incli- 
nation and the accretion disc truncation radius. We note, however, 
that these results depend stro ngly on the cont inuum spectral model 
and data reduction issues (see iNg et alfeOlOl for discussion). 

Although many aspects of AMP physics are known, one 
of the unresolved puzzles is th e origin of sudden pulse p rofile 
changes in SAX J1808.4-3658 JHartman et al] l2008l l200g) and 
other AMPs. There are many mechanisms that can cause such 
changes (Poutanen 2008) and they contribute to the "timing noise" 
seen in many AMPs (e.g.lBurderi et al.li2 006; Papitto et al. 2003; 
iRiggio et alj2008l ; |Patruno et al.l2009cl : |Poutanen et al..20 09). One 
of these remarkable pulse pro file changes occurred i n the 2008 out- 
burst of SAX J1808.4-3658 faartman et ai]|2009l) . In the begin- 
ning of the slow decay stage, the pulse profile was rather symmet- 
ric showing only small harmonic content. However, on Septem- 
ber 27 the fundamental pulse amplitude decreased by almost 50 
per cent, while no change was seen in other observed quantities, 
especially the observed flux and phases of the harmonics. A few 
days later on October 2 (MJD 54742), the pulse profile quickly 
morphed to have two peaks. This time the pulse profile transi- 
tion was accompanied by a decay in the observed flux and jumps 
in pulse phases. The double peaked profile was seen only for 
about three days, because on October 6 (MJD 54746), right be- 
fore the onset of the rapid drop stage, the fundamental amplitude 
increased significantly and the profile changed again to a shape 
similar to what was seen before September 27. This behaviour 
cannot be explained by a cha nge in the accretion dis c trunc ation 
radius that wa s sugg ested by Ibragimov & PoutanenI ( 120091) and 
IPoutanen et al.l ( 12009.) . Although this model explains the jump in 
the fundamental phase (and the pulse profile change), that is as- 
sociated with the rapid drop stage of the 2002 outburst of SAX 
J1808. 4-3658 (Burderi et al. 2006), it cannot be the case here sim- 
ply because the observed flux (which can be related to the trun- 
cation radius) in these sudden pulse profile changes is almost an 
order of magnitude higher a nd it remains nearly constant (see 
fig. 1 in Hartman et alj|2009l) . Among the 13 known AMPs (see 
IPatru no 2010a, and references therein), some sources show sim- 
ilar jumps in the fundamental pulse amplitudes as in the 2008 
outburst of SAX J1808.4-3658. A few examples include XTE 
J0929-314 dGallowav et alj2002h, the first eclip s ing AMP SWIFT 
J1749.4-2807 dMarkwardt & Strohmaverll2010l ; lAltamirano etal] 



I2OIII) and SWIFT J1756.9-2508 dPatruno et alj2010ah . Also, sim- 
ilar jumps in the fundamental pulse amplitude was seen in the 
2002 and 2005 outbu rsts of SAX J18 08.4-3658 (see fig. 6 in 
Ibragimov & PoutanenI, 20 09 and fig. 1 in lHartman et alj2008h . 

The variability in pulse amplitudes and profil es are some- 
times accompanied with jumps in pulse pha ses dPapitto et al.l 
I2OO7I ; iRiggio et al] |2008| ; IPatmno et al] l2009cb . In most cases, 
this timing noise is related wit h changes i n the observed flux 
dBurderi et al] |2006| ; llbragimov & Poutanei] |2009|; IPatruno et al] 
l2009d) and th e most notable changes usually occur in the funda- 
mental phase ('Burderi et al .11 20061 : [Papitto et alj|2007l ; IRiggio et alj 
[2008; Patruno et al. 2009c). This type of timing noise is most likely 
caused by mass accretion r ate dependent ch ange of the hotspot lo- 
cation on the NS surface iLamb et al J 1 2009.) , which explains the 
overa ll trends with flux and pulse phase residuals in SAX J1808.4- 
3658 (Hartman et al.l2008l) and in many other AMPs dPatruno et"al] 
|2009c). Howev er, alt hough the general trends can be explained 
dPatruno et alj|2009d) . there are several AMPs where additional 
timing noise is clearly present. Most notable cases are the 2002, 
2005 and 2008 outbursts of SAX J 1808.4-3658 XTE J 1807-294 
and XTE J0929- 314 (see fig. 3 in IPatruno et"ai]|2009d and fig. 1 
in H artman et al. l 2009). In these cases another mechanism - that is 
not related to changes in mass accretion rate - must be responsible 
for the timing noise. 

In this paper, we present our spectral and timing analysis of the 
2008 outburst of SAX J1808.4-3658 using Swift and RXTE data. 
We find evidence that the sudden changes in the pulse profiles and 
in pulse amplitudes of SAX J1808.4— 3658 are driven by a changing 
interaction between the NS magnetic field and the accretion disc. 



2 OBSERVATIONS AND DATA REDUCTION 

The 2008 outburst o f SAX J18 08.4-3658 was first detected on 
September 21 (Mark wardt & Sw ank 2008). The outburst was very 
similar to the 2005 outbur st in duration and bri ghtness, albeit the 
flaring stage was dimmer dHartman et al.ll2009l) . The slow decay 
stage of the outburst outburst was monitored by both Swift and 
RXTE from September 24 to October 3, during which the funda- 
men tal pulse amplitude a nd the pulse profile changed significantly 
(see iHartman et alj|2009l , fig. 1). In this paper we perform a de- 
tailed spectral and pulse profile analysis of these data in order to 
investigate the origin of this sudden timing transition. 

The RXTE data (ObsID 93027) were reduced using HEASOFT 
v.6.8 and the CALDB. We used data taken both by RXTE/PCA (3- 
25 keV) and HEXTE (25-100 keV). In cases where HEXTE expo- 
sures were short, we ignored the noisy channels above ~ 60 keV. 
Standard 0.5 per ce nt systematic was applied to the PCA spectra 
djahoda et alj|2006l) . To keep the calibration uniform, we used data 
from PCU 2 only (all layers). The Swift/XKT data (0.6-7 keV) were 
reduced using the XRTPIPELINE v.0.12.3 in HEASOFT v.6.8. The 
observations were performed in window-timing mode. We used 
standard filtering and screening criteria for the event selection. Ex- 
posure maps were generated with the XRTEXPOMAP task and the 
ancillary response files with XRTMKARF to account for different 
extraction regions (we used circular regions of 20 pixel radius), 
vignetting and psf corrections. The redistribution matrices (v.Oll) 
were taken from the CALDB. The XRT spectra were then grouped 
using GRPPHA to have at least 20 counts in each bin. 

We used the Swift/XKF data to model the time averaged spec- 
tra together with the quasi-simultaneous RXTE pointings. There 
were a few cases where we could not use all the XRT data. The 
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Table 1. The observation log. 



# 


ObsID 


Date and time 


Exposure" 


1 


93027-01-01-07'' 


2008-09-24 19:33:52-20:16:32 


2333/763 


2 


00325 827000^= 


2008-09-24 21:28:45-21:56:00 




3 


00030034026 


2008-09-24 23:17:40-23:32:00 


816 


4 


93027-01-01-06 


2008-09-25 05:24:48-06:24:32 


1035/414 


5 


93027-01-01-05 


2008-09-25 08:08:32-10:44:00 


6482/2359 


6 


00030034027 


2008-09-25 13:34:28-13:52:00 


1005 


7 


93027-01-01-10 


2008-09-25 19:52:16-20:36:32 


490/123 


8 


00030034028 


2008-09-26 11:59:49-12:19:00 


1142 


9 


93027-01-02-00 


2008-09-26 13:55:28-18:06:40 


9713/3264 


10 


93027-01-02-01 


2008-09-26 23:32:16-23:50:40 


1060/356 


11 


00030034029 


2008-09-27 04:04:15-04:23:00 


1116 


12 


93027-01-02-05 


2008-09-27 10:24:00-15:00:32 


8942/3411 


13 


93027-01-02-06 


2008-09-28 05:24:32-05:57:52 


1943/668 


14 


00030034030'' 


2008-09-28 12:11:35-15:44:00 


880 


15 


93027-01-02-03 


2008-09-28 16:13:52-17:35:44 


3201/1043 


16 


93027-01-02-04 


2008-09-29 06:22:08-07:45:36 


2775/1022 


17 


00030034031 


2008-09-29 10:44:25-10:59:00 


849 


18 


00030034032"= 


2008-09-30 01:12:05-01:33:00 




19 


93027-01-02-07 


2008-09-30 09:03:12-10:01:36 


3181/1154 


20 


93027-01-02-09 


2008-09-30 14:01:04-15:03:44 


2306/744 


21 


93027-01-02-02 


2008-10-01 11:44:48-14:59:12 


6637/2242 


22 


00030034033 


2008-10-01 14:25:42-14:45:00 


1131 


23 


93027-01-02-08 


2008-10-02 14:25:52-17:19:44 


6280/2035 


24 


00030034034-^ 


2008-10-02 22:23:50-23:59:59 


848 


25 


93027-01-03-00 


2008-10-03 09:49:20-13:22:40 


7557/2510 



°XRT or PCA/HEXTE exposure times (in seconds). ''Only data before 
X-ray burst used. '^Not used due to proximity to the X-ray burst. ''Only 
data from snapshot 2 used. "^Not used because the source was right on top 
of bad columns. -^Only data from snapshot 1 used. 

first XRT observation (observation 2, see Table [TJ was triggered 
by an X-ray burst. We found that the XRT spectrum differed from 
the PCA spectrum taken before the X-ray burst and also from the 
XRT spectrum taken only 3 hours after the X-ray burst (observa- 
tion 3). The spectra of observations 1 and 3 matched well indicat- 
ing that the first XRT spectra was affected by the X-ray burst and 
it was therefore excluded from the analysis. Also, we did not use 
the observation 18, because the source was right on top of the bad 
columns. In two occasions, where XRT observations were split into 
two snapshots (14 and 24), we used the data from the longer snap- 
shot. 



3 SPECTRAL AND TIMING ANALYSIS 

3.1 Simultaneous timing- and spectral transition 

During the slow decay stage of the 2008 outburst - on Septem- 
ber 27 (MJD 54736) - the pulse amplitude dropped a nd the pulse 
profile changed significantly (see lHartman et al.ll2009l fig. 1). The 
origin of these type of timing changes is not well understood, but 
they must be caused by some changes in the accretion geometry. In 
order to analyse this timing transition, we construct pulse profiles at 
various energies and different time intervals using ephemeris from 
iHartman et alj 120091). The observed pulse profiles are then fitted 
by a sum of two harmonics: 

F{<l)) = F{1 + ai cos[27r(0 - 0i)] + aa cos[47r(0 - ^a)]}, (1) 

where F((f)) is the flux at phase (j), F is the mean flux and ai, aa, (pi 
and (f)2 are the amplitudes and phases of the fundamental and the 



first overtone, respectively. The best-fitting amplitudes and phases 
are presented in Fig. [T] for the soft (3.7-5.7 keV) and hard (9.8- 
23.2 keV) energy bands. We selected these two bands because in 
the hard band we have emission only from the Comptonized com- 
ponent, whereas in the soft band there is a contribution from the 
blackbody component (e.g. Ilbragimov & Poutanenll2009f) . We see 
a drop in the fundamental amplitude and a change in the pulse pro- 
file on September 27 (MJD 54736), but we do not detect significant 
jumps in the pulse phases during the transition. 

However, we find that the timing transition was accompanied 
with a simultaneous softening of the energy spectrum below ~ 5 
keV. This can be seen by taking a ratio of the observed spectra be- 
fore and after the timing transition (from XRT observations 1 1 and 
14 and PCA observations 10 and 12). The ratio spectra are shown 
in Fig. |2] and the softening is clearly seen in the Swift as well as in 
the RXTEfPCA data. The fact that the timing- and spectral transi- 
tion occur simultaneously suggests a common physical origin, so 
a detailed broad band spectral analysis and pulse profile modelling 
of this transition is warranted. 



3.2 Spectral model 

We modelled the spectra using XSPEC 12 ( lAmaud|[l996l) . Errors 
are quoted at the 90% confidence level and the errors in the fluxes 
were computed with the CFLUX model in XSPEC. The reported lu- 
minosities are bolometric (calculated in the range of 0.01-500 keV 
from the best-fitting model) a ssuming a distance of D = 3.5 kpc 
JGallowav & Cummin j|2006l) . 

Our spectral model consist of CONSTANT x WABS x 
(DISKBB + BBODYRAD + DISKLINE + COMPPS). The spec- 
tral model is rather complex and has a large number of param- 
eters. The CONSTANT model component accounts for the differ- 
ent instrument normalizations between Swift/XKI, RXTEfPCA and 
HEXTE. We fixed the normalization of PCA and allowed XRT and 
HEXTE normalizations to vary. For XRT, the normalizations var- 
ied in a tight range around 0.9 between the different exposures, 
whereas HEXTE normalizations varied in the range of 0.52 — 0.62. 
The effect of interstell ar absorption was taken into a ccount us- 
ing the WABS model ('M orrison & McCammonI '1983Y which is 
parametrized by the absorption column A'h- We m odelled the ac- 
cretio n disc component with the DISKBB model JMitsuda et al.l 
119841) . which has two parameters; the inner disc temperature Tdisc 
and normalization (proportional to inner disc radius i?disc) ^disc ~ 
[(-Rdisc/kin)/(Z)/10 kpc)]^ cosi. For the thermal emission from 
the NS surface, we used the BBODYRAD model, which is charac- 
terized by the black body temperature Tbb and normalization (pro- 
portional to the surface area) Khb = [(-Rbb/km)/(Z)/10 kpc)]^, 
where i?bb is the appa rent black body radius as obs erved from in- 
finity. Model COMPPS jPoutanen & Svenssorjl 19961) - characteriz- 
ing Comptonization in the accretion column (we assumed slab ge- 
ometry) - is described by the following parameters: Thomson op- 
tical depth r across the slab, electron temperature Tc, seed photon 
temperature Tsccd for Comptonization, and normalization (propor- 
tional to the surface area) Ksacd ~ [(-Rscod/km)/(-D/10 kpc)]^, 
with -Rscod corresponding to the size of the accretion column. The 
Compton reflection component (included in COMPPS) is defined by 
the amplitude 5R = Q./2-n (where Q. is the solid angle covered by 
the cold reflector as viewed from the X-ray emitting hotspots), in- 
dex a of the radial profile of the disc emissivity (oc r"), inner- and 
outer disc radius, inclinatio n i and the iron lin e energy and nor- 
malization iCpc (DISKLINE. iFabian et"al]|l989l) . We note that the 
DISKLINE model inner disc radius can be substantially different 
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Figure 1. Change of the pulse profile parameters during the outburst. Left panels are for the soft (3.7-5.7 keV) band, while the right panels are for the hard 
(9.8-23.2 keV) band. Upper panels present amplitudes of the fundamental (red squares) and of the first overtone (blue diamonds). Lower panels correspond 
to the pulse phases shown by the same symbols. Only those phases and amplitudes that can be constrained are plotted, and errors for only a few representative 
data points are plotted for visual clarity. The time of the spectral transition is marked with the grey shaded strip (see Fig.[2j. Note the absence of a pulse phase 
shift during the transirion (MJD 54736). 
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Figure 2. Ratio of the spectra from observations after and before the timing 
transition. The coiresponding XRT spectra are from observations 14 and 
11 (blue squares) and the PCA spectra are from observations 12 and 10 
(red triangles). The spectral data are binned for visual clarity. The spectral 
transition occurred between 2008-09-27 04:23:00 and 2008-09-27 10:24:00 
(see Table[T). The increase of the flux is most evident below ~ 3 ke V, where 
the accretion disc component dominates (see lPatruno et al.ll2009al and Fig. 
0. 



from i?disc, which we derive from the DISKBB model normaliza- 
tion. The former is related to the disc truncation radius Rt, whereas 
the latter traces the radius where the dissipation ceases, which is re- 
lated to the magnetospheric radius i?m (see end of Section [331 >. 
When we attempted to fit all these parameters to the data, 



we found several unphysical changes in the best-fitting parameters 
throughout the outburst. Most notable inconsistencies were variable 
inclination i and a larger hotspot radius than a typical NS radius. 
Furthermore, several parameters showed tight correlated changes 
among them. Tc, r and 5R - which together define the spectral slope 
above ~10 keV - were strongly correlated, resulting in large uncer- 
tainties in these parameters. Also, A^h showed inconsistent changes 
during the outburst. We thus had to fix several parameters to obtain 
reasonable values. We did this in an iterative manner, where we 
fixed certain parameters, then fitted the data again, inspected the 
results and then fixed more parameters if necessary. 

Initially we fixed five parameters. Incl ination was fixed toi = 
60° ( for both reflection and iron line, as in llbragimov & PoutanenI 
|2009). The inner disc radius that affects both the smearing of the 
line and Compton reflection was fixed to a value of 10i?s, where 
Rs = 2GM/(? is the Schwarzschild radius. The outer disc radius 
was fixed to 1000i?s- The radial disc emissivity for reflection was 
fixed to Of = —3. Also, because of poor energy resolution of PCA 
and poor statistics of XRT spectra in the iron line region, we fixed 
the iron line energy to 6.4 keV. 

The model as such described the data very well (the best- 
fitting parameters and fitting statistics were compatible to the 
"final results" present e d in Table [2^. However, as discussed i n 
iGierliiiski & PoutanenI ( |2005|) and llbragimov & PoutanenI ( |2009|) , 
the hotspot surface area (ex Khh) was poorly constrained and 
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in many cases resulted in large and unphysical values. Further- 
more, some residuals were present in the XRT data around 2 keV, 
where the calibration is known to have problems because of the 
Si Ka edge at 1.839 keV (Godet et al. 2009). These residuals cre- 
ated a problem in the fitting, because the flux contribution of the 
Tbb ~0.5-0.6 keV NS blackbody component and the Tdisc ^0.3- 
0.4 keV disc component was very similar at this energy. We identi- 
fied two cases (6 and 11) - where the XRT observations were about 
3' off-axis - in which these instrumental residuals were the main 
reason why Tdisc, ^disc, Thh, Khb and A^h seemed to change in 
a correlated way during the outburst. The best-fitting value of A'^h 
varied roughly in a range of (0.9-0.95)±0. 15 x 10^^ cm~'^ for most 
of the observations, except in these two cases where it dropped to 
a value of ~ (0.7 ± 0.05) x 10^^ cm~^. As we do not expect A^'h 
to vary and its value was tightly correlated with the parameters of 
the accretion disc (Tdisc and ii'disc), we fixed this parameter to the 
Galactic value of 1.13 x 10^^ cm"^ (obtained from the HEADAS 
tool NH) to avoid spurious chang es in the disc parameters. Also, 
following llbragimov & PoutanenI (12002), we set Khb = J^^scod in 
the fitting, which helped to remove the correlation between Tbb and 
Tdiac described above. 

After another round of spectral fitting, we found that 
Tc ~(45-55)±10 keV in all the spectral fits. In some cases the 
spectra were very noisy above ~60 keV and some small changes 
in Te and r between different observations were identified to be 
caused by one or two spectral channels only. As we did not see any 
statistically important trend between the observations, we fixed Te= 
50 keV, which was the best-fitting constant value. We then fitted the 
data again, and got physically acceptable values for all the model 
parameters. We also checked, that the removal of XRT data in the 
1.6-2.4 keV range did not change the best-fitting parameters sig- 
nificantly. 



3.3 Results of the spectral modelling 

The time evolution of the best-fitting parameters of the spec- 
tral modelling are presented in Fig. [3] and Table |2] It is imme- 
diately obvious that there is a clear transition in the best-fitting 
spectral parameters on September 27 (MJD 54736), exactly when 
the timing transition also occurs. However, we emphasise that 
the bolometric luminosity during this transition does not change 
(Lhoi ~ 8 X 10^® ergs~^ see Fig.jS] top panel). The most promi- 
nent changes are in the parameters of the accretion disc. In the 
transition, the disc component becomes more luminous, the inner 
disc temperature Tdisc increases and the inner disc radius i?disc de- 
creases. Before the transition, the bolometric disc luminosity was 
Ldisc ~ 2.4 X 10'^^ ergs~^ and after the transition it increased 
to Ldisc ~ 3.2 X lO'^^ergs^^. This means that the luminosity 
of the accretion disc changed from roughly 30 to 40 per cent of 
the bolometric luminosity, whereas the luminosity of the hotspot 
decreased. We also detect a significant increase of the hotspot tem- 
perature (and the seed photon temperature) and a decrease in the 
apparent radius of the hotspot. We also note that we do not see any 
change in the optical depth of the shock as r ~ 1 throughout the 
outburst (see Table|2}. This suggests that although the characteristic 
temperatures and radii change in the transition, the properties in the 
Comptonized emission above ~ 10 keV remain constant. We make 
use of this result in Section |4~2l where we analyse the simultaneous 
pulse profile transition. 

The spectral fitting indicates that -Rdisc decreases from ~50 
to «30 km in the transition, while the inner disc temperature Tdisc 
increases from ~0.3 to ~0.4 keV. This effect is seen in both Swift 
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Figure 3. Time evolution of the spectral parameters. Luminosities are given 
in units of 10'"' ergs^^, temperatures in keV and radii in km. The error 
bars in time axis denote the time between Swift and RXTE pointings. The 
grey ai'ea denotes the time period when the spectral and timing transition 
must have occurred (see the text and Fig.[2)- The spectral transition stands 
out as an increase of the inner disc temperature Tdisc, a decrease of -Rdisc 
and a change in the hotspot radius and temperature. The bolometric lumi- 
nosity remains roughly constant at Lbol ~ 8 X 10'^'' ergs"^, which in- 
dicates a constant mass accretion rate M through the transition. Note also 
that the values of -Rdisc and Rhh are affected by systematic uncertainties, 
see text. 



and RXTEfPCA spectra, which is illustrated in Figures|2l[3]and|4] 
Interestingly, after the transition we see that -Rdisc ~ Rco, as the 
co-rotation radius 

7?co = (GA/./ajJ)'/^ (2) 

of SAX J1808.4-3658 is about 31 km for a NS mass of M. = 
1.4Mq (uj, is the angular spin frequency of the NS). There is of 
course an inherent uncertainty in deriving -Rdisc as its value de- 
pends on the assumed distance and the inclination as -Rdisc oc 
D/\/cosi. However, the distance is known rather accurately 
JGallowav & Cumminel 20061) and we will show in Section |4!2l that 
our initial assumption of i = 60° is consistent with the pulse pro- 
file modelling. Other complications in relating modelled radii (both 
-Rdisc and -Rbb) to "physical" radii are related to the deviation of the 
local spectra from the blackbody, assumed boundary conditions and 
relativistic effects. Because these effects do not change the interpre- 
tation of the results qualitatively - in the sense that there is only a 
systematic shift in the modelled radii - we defer the discussion of 
these effects to SectionlSTSl 
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Table 2. Best-fitting parameters. The fixed parameters were: JVjj = 1.13 X 10 



:50keV, i = 60°. 



#" 



-^ disc 

(keV) 



-'^-disc 

(km/10kpc)2 



Tbb 

(keV) 



fxlO" 



-^ seed 

(keV) 



(km/10 kpc)2 



-'^ disc 



380t??» 
430tigS 
410t^go 
450t?^S 

180l«° 
280t^S° 
SOOt^go 



xL/d-o.f. 



1,3 


0.28 


5,6 


0.29 


8,9 


0.28 


10, 11 


0.29 


13,14 


0.40 


16, 17 


0.42 


21,22 


0.33 


24,25 


0.32 



+0.02 
-0.02 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.02 
+0.01 
-0.02 
+0.02 
-0.02 
+0.01 
-0.01 
+0.01 
-0.01 



1250l^?» 
1070111° 
llSOt^^o 

40011° 
320lig 
650l°° 
59011J0 



0.53 
0.51 
0.52 
0.52 
0.64 
0.66 
0.55 
0.52 



0.07 

0.10 
+0.06 

0.06 
+0.04 

0.05 
+0.07 

0.12 
+0.09 

0.09 
+0.07 

0.08 
+0.06 

0.06 
+0.06 

0.07 



1.50 
1.50 
1.55 
1.62 
1.51 
1.59 
1.44; 
1.30 



+0.36 

0.35 
+0.25 

0.25 
+0.23 

0.22 
+0.44 

0.43 
+0.46 

0.45 
+0.40 

0.40 
+0.22 

0.22 
+0.18 

0.18 



0, 

0.84 

0.84 

0.82 

1.00 

i.oa 

0.87 
0.80 



0.07 

0.09 
+0.06 

0.06 
+0.04 

0.05 
+0.07 

0.11 
+0.09 

0.09 
+0.08 

0.09 
+0.06 

0.06 
+0.06 

0.07 



1.05 
1.05 
1.05 
1.09 
1.05 
1.09 
1.02 
0.98 



+ 0.11 
0.08 

+0.06 
0.05 
0.05 
0.05 

+0.08 
0.11 
0.08 
0.08 

+0.12 
0.09 
0.08 
0.06 

+0.10 
0.08 



< 0.40 



+0.11 
0.12 

+0.11 
0.10 



0.21 
0.21 

< 0.34 

< 0.42 

< 0.35 

27+°-i^ 
'-'■^'-0.16 

0.38toj? 



1.60 
1.70 

1.47 

l.i 



2.17 
2.17 
1.68 
1.29 



+0.08 

0.07 
+0.05 

0.05 
+0.04 

0.04 
+0.05 
-0.06 
+0.07 

0.08 
+0.08 

0.08 
+0.05 

0.05 
+0.04 

0.04 



5.15 
5.37 
5.04 
5.29 
5.53 
5.38 
4.42 
3.52 



+0.17 

0.14 
+0.10 

0.09 
+0.08 

0.08 
+0.15 

0.17 
+0.25 

0.14 
+0.19 

0.15 
+0.10 

0.08 
+0.09 

0.07 



1.07/371 

1.21/550= 

1.07/545 

1.22/ 553^^ 

1.01/535 

1.08/533 

1.11/533 

1.02/484 



" Observation numbers from Table [T] * Flux in units of 10 " erg cm 



■ Poor fits are mainly due to large residuals in XRT spectra around Ri 2 keV, 



which are most likely due to a calibration issue at the Si Kq edge at 1.839 keV JGodet et aL ,2009i). Removing the data between 1.6-2.4 keV improves the fits 



tox^c 



1.1, without a significant change in the best-fitting parameters. 



In accretion on to neutron stars one expects that a constant 
fraction of the available accretion power is released in photon lu- 
minosity; what is not radiated in the disc, will be radiated at the 
NS surface. Therefore, the observed luminosity is thought to be a 
good measure of the mass accretion rate M . As the bolometric lu- 
minosity during this transition does not change, we can conclude 
that the sudden change in the disc parameters cannot be caused by 
a sudden increase in the mass accretion rate. This points towards 
the conclusions that the change in i?disc is not related to a change 
in the truncation radius (which is proportional to the Alfven radius, 
where the pressure of the NS magnetic fiel d equals the ram pressure 
of the accreting gas, e.g. lLong et ahluOOSl) 



Rt =kAi2GM,y 



-2/7 4/7 



(3) 



Here /i — BsR'^ is the NS magnetic moment and 7?* is the NS 
radius. The truncation radius Rt is also commonly called the mag- 
netospheric radius, but in the following we instead use this term 
to define the outermost radius Rm where the field lines of the NS 
are coupled with the accretion disc ( Lovelace et al. 1995, hereafter 
LRB95). We propose that the interaction region between the disc 
and the NS magnetic field is large (so that Rm > Rt, see Fig. 
131 bottom panels) and that the modelled inner disc radius actually 
corresponds to the magnetospheric radius /Jdisc ~ ^m and not the 
truncation radius Rt- 



4 VARYING DISC-MAGNETOSPHERE COUPLING IN 
SAX J1808.4-3658 

4.1 Physical picture 

We interpret the spectral transition to be caused by a change in the 
way the dipole field of the NS is coupled with the accretion disc. 
The reason for the observed spectral and timing transition is most 
likely caused by opening of the field lines. 

Many theoretical stud ies of the disc-magneto sp here coupling 
have been published (e.g. iGhosh & Lamb 1979a b; Wang 1 9871: 
Lovelace et alJl995l : lRapDaDort et al.2004. ; .Kluzniak & RappaportI 



20071). Here we use the framewo rk developed by LRB95 to inter- 
pret the observed transition (see Uzdenskvl l2004l for a review of 
alternative models). In the LRB95 model, the magnetic field con- 
figuration is such that close to the accreting NS the field lines are 
closed (defined here as the magnetosphere) and outside the mag- 
netosphere the field lines are open (see Fig. |4ll. The field lines de- 
couple from the disc because of a large difference in the angular 



velocity of the Keplerian disc and the NS (except at Rco)- As the 
disc is a good conductor, the NS magnetic field is "frozen" into it 
and the twisting of the field lines caused by the differential rota- 
tion creates a toroidal field component B^ out of the poloidal field 
Bz- However, there is a critical twist of 7c — B^/Bz ^ 1 after 
which the increased magnetic pressure tends to open the field lines, 
thus severing the link b etween the NS a nd the disc (see LRB95, 
lUzdenskv et alj2002al lbl: IUzdenskvll2004l and references therein). 

Within the magnetosphere (r < Rm), the magnetic field 
threading the disc can force the gas to co-rotate with the NS so 
that j jjjr) ~ a;« (LRB95), w hich is also seen in MHD simulations 
(e.g. iRomanova et alj|2002h . This has a strong effect on the radia- 
tion flux emitted from the disc, which is most easily seen from the 
energy conservation equation (see eq. 7 in LRB95) 



S/^t ( r— ) + — 
dr 



dzfjtJ = 2crsBT'cff. 



(4) 



Here E is the surface density, Vt is the turbulent viscosity (which 
is assumed to be describe d by the alpha prescription Ut = ac^H, 
IShakura & Sunvaevll 19731) , Cs is the sound speed, z is the coordi- 
nate along the normal to the disc plane, H is the disc height, rjt is 
the magnetic diffusivity of the disc, J is the current density, ctsb is 
the Stefan-Boltzmann constant and Tcff is the effective temperature 
of the disc. In equation (|4j the first term in the left hand side is the 
viscous dissipation rate per unit area, the second term is the Ohmic 
dissipation rate and the right hand side is the flux radiated by the 
optically thick and geometrically thin disc. If Lj(r) ~ cj, within the 
magnetosphere (r < i?m), the viscous dissipation ex (dcj/dr)^ is 
strongly suppressed and the Ohmic dissipation becomes the dom- 
inant energy release mechanism (LRB95). Expressing the current 
density through the accretion velocity Vr and vertical magnetic field 
strength B^ = Bs{R*/r)^ as J = {c^ /4n'r]t){vr/c)Bz and inte- 
grating equation Q over the radius we get an order of magnitude 
estimate of the lu minosity due to the Ohmic dissipation (see also 
IWangetalJI 19901) : 



Loh 



H 



Rt 



Hvj. 
Vt 



^ *Rt 



(5) 



For geometrically thin disc H/r ~ 0.01 with Ca < 10* cm s~^, 
substit uting Vr ^ aCsHIr (e.g . iFrank et al.ll2002h . rit '^ ft ~ 
aCsH I IFromang & Stond 120091) and a < 1, we can show that 
iohm < 10^^ ergs~^ forthe typical parameters of SAXJ1808.4- 
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Figure 4. Photon spectra before the transition (from observations 10 and 11, see Table[T) and after the transition (from observations 13 and 14) are shown in 
the top panels (a) and (b), respectively. The blue, red and green data points are for Swi/f/XRT, RXTE/PCA and RXTE/UEXTE, respectively. The solid black 
line shows the model spectrum. The pink dot-dashed Hne is the unabsorbed disc component and the three-dot-dashed line the absorbed disc component. The 
purple dotted line is the NS blackbody component, whereas the dashed line is the Comptonized component. The weak Compton reflection and the iron line 
are below the plotting range so that the spectral transition is more easily visualized. Pulse profile before the transition (from observations 9 and 10) and after 
the transition (from observation 12) are shown in the middle panels (c) and (d), respectively. The solid Une is the contribution from both spots, whereas the 
dotted line shows the contribution of the "main" spot. Lower panels illustrate the geometry of the disc-magnetosphere interaction region before (panel e) and 
after (panel f) the transition. The dark region in the accretion disc within Rt < r < R^ is assumed to rotate as id{r) fs w, so that the energy dissipation rate 
in that region is strongly suppressed. 
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3658 (Bs ~ 10* G, R, ^ 10 km, Rt ^ 20 km)QThus the Ohmic 
dissipation within the magnetosphere is much smaller than the vis- 
cous dissipation outside the magnetosphere (r > Rm), and there- 
fore does not contribute much to the observed flux from the disc. 
Therefore, the inner disc radius -Rdisc derived from the spectral fits 
actually should correspond to the magnetospheric radius Rm- This 
makes the physical interpretation of the transition rather straight- 
forward. 

If -Rdisc ~ ^m, the reconfiguration during the transition can 
be understood as opening of the field lines (see Fig. [4] bottom pan- 
els). Initially the interaction region between the magnetosphere and 
the disc is large, such that it extends from the truncation radius 
Rt ~ 20 km to Rm ^ 50 km. Then the field lines open resulting 
in a smaller Rm ~ 30 km, so that the region of the suppressed 
dissipation (where a; ~ oj*) is smaller (dark regions in the disc in 
the bottom panels of Fig. 0. This readily explains why there is an 
increase in the inner disc temperature Tdisc, a decrease of -Rdisc and 
why -Ldisc/iboi increases. 

The fitting indicates that -Rm ^ Rco before the transition, but 
we note that the exact value of Rm is rather uncertain (see Sec- 
tion |5.3l ). According to Rappaport et al.l <2004) Rm can exceed the 
corotation radius by 30 per cent before the centrifugal barrier halts 
the accretion. This might set the upper limit for -Rm before the 
transition. Then, the sudden opening of the field lines (decrease 
of -Rm) could be caused by a change in the properties of the ac- 
cretion disc. In the LRB95 model -Rm is determined by the critical 
twist 7c = BJB^ = -Hr{uj{r) - uj*)/rjf As 7c « 1 (e.g. 
lUzdenskv et alj|2002al lbl) is thought to be a constant, a change in 
the magnetospheric radius Rm - for a fixed M - could be caused 
by a change in the magnetic diffusivity rjt . Furthermore, as the dif- 
fusivity r]t is most likely caused by the sa me MHD turbulence tha t 
produces the viscosity Ut (so that -qt ~ i^t. lFromang & Stonel2009r) 
we speculate that the observed spectral- and timing transition could 
be related to changes of viscosity in the accretion disc. Also, re- 
cent MHD simulations have shown that if the magnetic Reynolds 
number (= CsH/rit) is below a certain critica l value, the magneto- 
rotational instability l IBalbus & Hawlevll998h that is believed to be 
responsible for the angular momentum transport in the accretion 
disc might quench, possibly causing high- and low viscosity states 
([Simon et al. 201 1). If the pulse profile variability - and in general 
the timing noise - in AMPs is ultimately caused by such a mech- 
anism, observations of AMPs could be used to place constraints 
for rjt and Vt because for many AMPs the key parameters such as 
-Bs , M and the relevant radii and temperatures are rather accurately 
known. 



4.2 Pulse profile modelling 

The observed spectral transition can be caused by a change in 
the coupling between the magnetic field of the NS and the ac- 
cretion disc. Changes in the geometry of the disc-magnetosphere 
interaction region, on the other hand, alter the way the accretion 
flow gets channelled onto the NS surface and cause variations in 
the hotspot size and shape, which affect the pulse profiles (see 
Fig. El). MHD simulations of i Romanova et alj J2004I ) showed that 
the spot shape strongly depends on the magnetic inclination 6, 
the angle between the magnetic pole and the rotational axis. For 
e < 15° - which is most likely the case for SAX J1808.4-3658 



^ We constrain the truncation radius to -Rt ~ 20 km from the pulse profile 
modelling in Section l431 




Figure 5. Geometry of the hotspots. The magnetic dipole is inclined by 
angle 8 to the spin axis and the co-aligned hotspots are assumed to be ring- 
shaped. The inner edges of the hotspots are displaced from the magnetic 
poles by angle pm and the outer edges by pout ■ The accretion disc trunca- 
tion radius, that determines which part of the secondary hotspot is blocked, 
is related to pout using equation (6). The pulse profile vaiiations can be 
modelled simply by changing Pm (which is a function of -Rm) as shown in 
SectionlO 



(e.g. Ilbragimov & Poutanenl 120091) - the spot shape is close to a 
ring. 

We thus assume that the hotspot shape is a ring around the 
magnetic pole so that the size of the hotspot is determined by the 
inner hotspot edge pm and the outer hotspot edge pout (see Fig. 
[5]l. By further assuming that the NS magnetic field is a dipole, pin 
and Pout can be related t o the magnetospheric radius and the dis c 
truncation radius (see e.g. lFrank et al.ll2002l ; |Poutanen et alj|2009n : 



(6) 



sin Pin = cos9\/Rt/Rm, sinpout = cos9^R,/Rt. 

We also assume that changes in the magnetosphere size do not af- 
fect the truncation radius -Rt . This can be justified as the detailed 
mechanism of the interaction should not play a big role owing to 
the strong (oc r"^) radial dependence of the magnetic energy den- 
sity, which ~ for a given M - determines the location of -Rt . Ob- 
viously parameters such as Al,, R,, i, and -Bs do not change in 
the transition either. What is expected, however, is that the hotspot 
size should change, which was also seen in the modelling of the 
time averaged spectra in section [373] This can be seen as a natural 
consequence of the opening of the field lines. As illustrated in Fig. 
[4] the field lines that connect the disc at -Rm are the ones that con- 
nect closest to the magnetic pole of the NS. Given our assumption 
that Rt is constant, pout is also constant as is constant. Therefore, 
under these assumptions, the change in the hotspot size can only be 
caused by a change of the inner hotspot edge pin (see Figs |4] and 
[5J. Furthermore, the opening of the field lines would not change 
the hotspot centroid, and thus would not cause a jump in the pulse 
phases. As we do not see any such change in phases immediately 
after the pulse profile change (see Fig.[T), our choice of a hotspot as 
a symmetric ring around the magnetic pole seems to approximate 
the emission region to a sufficient accuracy. 

We use the PC A data from observations 9 and 10 for the 
"pre-transition" pulse profile and observation 12 for the "after- 
transition" pulse profile. We only concentrate on the pulse profile 
changes in the hard band (9.8-23.2 keV) for two reasons. (1) In 
the hard band we have emission only from the Comptonized com- 
ponent, whereas in the soft band blackbody component also con- 
tributes. (2) In the hard band we have emission from high scatter- 
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ing orders jViironen & Poutanenll2004l : |Poutanenll2008h . where the 
energy and angular distributions of the radiation intensity I{E, () 
at the stellar surface can be approximated with a simple formula: 



Table 3. Best-fitting parameters from the pulse profile modelling. 



I{E,() oc {l-hcosC)E' 



-(r-i) 



(7) 



where E is the photon energy, (" is the angle relative to the sur- 
face normal, h is the anisotropy parameter and F is the power-law 
photon index (T ^ 2 for SAX J1808.4-3658). 

We fitted the profiles with the model describe d in de- 
tail inlPoutanen & G ierlihski (2003) and Poutanen & Bel oborodovl 
( 12006 "). The model accounts for special relativistic effects (such 
as Doppler boosting and aberration) as well as the general rela- 
tivistic effects (redshift and light bending). We also account for 
the eclipses of the hotspot by the accretion disc as described by 
llbragimov & P outanen (2009). We initially arbitrarily chose Af, = 
1.4M0, R, = 10.3 km, i — 60°. We set 9 and pout as free param- 
eters, but forced them to be the same for the two pulse profiles. Rt 
was computed from pout and 9 using equation (|6}. An additional 
free parameter is the phase shift <jf>shift describing the displacement 
of the pulse peak. Because the extracted pulse profiles were sep- 
arated by one day and no shifts in phases were seen (see Fig. \Q, 
we force (j>BhUt to be the same for both pulse profiles. The only 
parameters that we initially let to vary in the transition were pin 
and h, which in our view could change if the disc-magnetosphere 
coupling changes. 

However, the strong secondary pulse after the transition 
(which we associate with the secondary spot) proved to be prob- 
lematic to model. The only way we could get such a peculiar shape 
was that the secondary spot was partially blocked by the accretion 
disc (we assume that the secondary spot is antipodal). We realised 
that the reason we did not get a good fit to the data was our guess of 
M«, R, and i. We then let R, and i to vary in the fitting and found 
that there was a rather tight range of these parameters that we could 
produce such a pulse shape. We found that for M* = 1.4M0, only 
with a NS radius of _R* ~ 11 km were we able to produce the pulse 
shape correctly. Furthermore, in the initial fitting we did not see a 
significant change in the h parameter between the two pulses as it 
varied in a tight range of /i ~ 0.57±0.05. This is actuall y expected 
since h fundamentally depends on the optical depth r jPoutanerj 
2003), which did not change in the transition (see Fig. [3] and Table 



2)1. Va lues of /i > 0.5 correspond to optically thin shocks JPoutanerj 
20081) ■ so the values of r ~ 1 and h ~ 0.57, obtained here indepen- 



dently from spectral and pulse profile modelling are consistent with 
each other. Therefore, we let h to be free parameter, but required it 
to be the same for both pulses. 

Thus, the modelling was performed with parameters h, pout , i, 
9 and <j)sMit and the parameters that were allowed to vary between 
the pulses were pin,i and pin, 2- The best-fitting parameters of the 
modelling are presented in Table [3] The corresponding pulse pro- 
files are shown in Figs|4jc) and (d), where the dashed line shows 
the contribution of the main hotspot, and the solid line shows the 
sum of the main and the secondary hotspots. The results indicate 
that the secondary hotspot is almost entirely blocked (by the NS 
itself and the accretion disc) around phase 0, and becomes partially 
visible around phase 0.5. The overall contribution of the secondary 
hotspot to the observed flux increases in the transition from about 
5 to 10 per cent mostly because of the increase of the hotspot tem- 
perature (associated with the decreasing size). We also see that for 
the assumed M* = 1.4M0 and _R, = 11 km, we get physically 
plausible results for all the free parameters. 

(i) The inclination is roughly i ~ 60°, which is consistent with 



Parameter 


Value 


Units 


NS mass M* 


1.4 (fixed) 


Me 


NS radius i?* 


1 1 (fixed) 


km 


Inclinadon i 


58t^ 


deg 


Magnetic inclination 9 


11±1 


deg 


Outer spot radius pout 


46 ±6 


deg 


Inner spot radius before transition Pin,i 


<12 


deg 


Inner spot radius after transition pin.2 


31 ±2 


deg 


Disc truncation radius Rt '^ 


211^ 


km 


Anisotropy parameter h 


0.57 ±0.03 




X^/d.o.f. 


62.7/57 





Rt is not a free parameter, it was computed using Eq. (6). 



previ ous estimations jDelove et al.l 120081 : llbragimov & PoutanenI 
I2OO9I1 . 

(ii) The value of the anisotropy parameter h « 0.57 is also what 
we ex pect from Comptonization of optical depth t « 1 iPoutanenl 
l2008h . that was seen in the time averaged spectra. 

(iii) The magnetic inclination 9 « 10° is consistent with 
the a ssumption that the s hape of the hotspots are close to a 
ring I Romanova et al. 2004) and similar to previous constraints 
( llbragimov & Poutanen.2009.) . 

(iv) The truncation radius Rt ~ 20 km, which is within the co- 
rotation radius of R^o = 31 km and is smaller than R^ ~ R^o, 
which is what we expected. 

(v) Most importantly, the increase of pin from < 12° to ~ 30° 
is what is expected, if the transition in the spectrum and the pulse 
profile are both caused by a decrease of Rtn- Therefore, we see that 
the results of the pulse profile modelling are consistent with our 
interpretation and - taken together with the spectral information 
- fully support the picture that the pulse profile changes in SAX 
J 1808.4-3658 are a consequence of a dynamic magnetosphere. 



5 DISCUSSION 

5.1 Implications for spin-up torques during the outburst 

The changes in the interaction between the disc and the magne- 
tosphere can have a large effect on the torque that spins up (or 
spins down) the NS. Several different models have b e en developed 
to compute the torques (e.g.'Ghosh & Lamb^'l979a''D;'Wans" 19871; 
Lovelace et al. 1995; Rappaport et al. 2004; Matt & Pudritz 20ol; 
iKluzniak & RapDaDortll2007l : lMatt et al.ll2010l) . General ly, in addi- 
tion t o the "accretion torque" Ta = M^GM^Rt (e.g. iMatt et al.l 
I201CI) . the disc-magnetosphere coupling can cause an additional 
"magnetic torque," which is strongly dependent on the nature of the 
interaction and varies significantly between the models. Neverthe- 
less, it can be concluded that small variations of Rm independent of 
the mass accretion rate (as observed here for SAX J1808. 4-3658) 
can cause the NS spin-up or spin-down depending on the geometry 
and position of Rm relative to _Rco- 

Also the data suggests, that only those models that predict that 
oj ~ oj* in the magnetosphere can be valid to explain the observed 
spectral transition. So observations of AMPs - as presented here - 
can be used to test different hypothesis and constrain the nature of 
the disc-magnetosphere coupling. 
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5.2 Implications for the timing noise in AMPs 

Although the pulse amplitude dropped significantly, the pulse phase 
remained constant during the transition. This is actually not sur- 
prising, based on the fact that the magnetic inclination seems to 
be about 6 ~ 10°. As the ho tspot shape is most likely close 
to a ring JRomanova et al.ll2004h . the change in _Rm should only 
affect the location of the inner edge of the hotspot pin, but not 
change the location of the hotspot on the NS surface. How- 
ever, the situation might be different for other AMPs, where 9 
can be larger. For su ch AMPs, the hotspot shapes are most cer- 
tainly not symmetric JRomanova et alj|2004h . and changes in Rui 
might result in a change in the shape, size, latitude and -most 
importantly- the hotspot longitude. Shifts in the hotspot longi- 
tude can especially cau se large jumps in the observed pulse phases 
(see iLamb et alj |2009| , for discussion). We speculate that such 
Jim-dependent motion of the hotspots (and the associated phase 
jumps) could be the origin of at least part of the timing noise in 
AMPs. Observationally, such mechanism could cause the outliers 
in the X-ray flux - phas e residual relations i n XTE J 1807-294 



Riggio et al.|[2qo3:|Patruno et al...20093.l2010b) . XTE J0929-314 



Galloway et al.l2002l:|Patruno et al.l2009d) and IGR J 1751 1-3057 
Riggio et al.ll201 lUlbragimov et alj|201 ih . whereas the overall X- 



ray flux - phase residual trends can be caused by M-dependent 
motion of the hotspots JLamb et all l2009l : iPatruno et alTl2009d. 
l2010br ). Also, we stress that an additional contribution to the tim- 
ing noise can be caused by A/-dependent variation of the trun- 
cation radius R t, as it affects the visibility of the secondary spot 
jlbragimov & P outanen 2009; Poutanen et al. 2009). These factors 
do not however exclude the possibility that some AMPs do spin- 
up during X-ray outbursts (as is most likely the case with IGR 
J00291+5934. e.g.lFalanga et al.l20oilBurderi et al.l2007l : |Patrimol 



1201 Obi : iHartman et al J201 ll : |Papitto et al J201 ih . 



5.3 Possible caveats 

(i) The truncation radius of Rt ~ 20 km was obtained from 
pulse profile modelling. We assumed that this radius is not affected 
by the change in the disc-magnetosphere coupling that we propose 
to be the reason for the transition. We also assumed that the outer 
edge of the ring-shaped hotspot pout is determined from equation 
(|6j, which is not necessarily the case as the field topology should 
differ from a pure dipole. These factors can change the constraints 
on the other parameters. However, a detailed anal ysis of these fac- 
tors s hould be done by 3D MHD simulations (as in lRomanova et al.l 
[20041), but such a study is clearly beyond the scope of the present 
paper. 

(ii) As the Ohmic dissipation in the magnetosphere is small, we 
proposed that -Rdisc ~ Rm- However, there are some uncertainties 
in deriving -Rdisc in addition to those related to spectral modelling 
in Section [T2I The main systematic uncertainty in estimating Rm 
from -Rdisc comes from the fact that Tdisc is in reality a colour 
temperature and not the effective temperature of the disc. Also, the 
inner disc boundary condition might differ significantly from what 
is assumed in the DISKBB model, where the radial dependence of 



temperature is T{r) 



<i/i 



These effects could be, in principle. 



accounted for to obtain the corrected inner disc radius 



-^disc sjcol-^disi 



(8) 



where ^ accounts for the fact that the modelled inner disc temper- 
ature Tdisc does not actually occur at - Rdisc and /coi is the colour 
correction factor JGierlinski et al J 19991 ). The colour correction fac- 



tor was computed to be /coi ~ 1.7 bv lShimura & Takaharal j 19950 
for black hole systems in the soft state and it seems to be w eakly de- 
pendent of the mass accretion rate (see Davis et al. 2005). T he cor- 
rection factor ^ = 0.37 was computed by Gierlihski et al. j 19991) 
by assuming a zero torque boundary condition at the inner edge of 
the disc in the pseudo-Newtonian potential. However, because of 
the different inner disc boundary conditions for accretion onto an 
AMP 5 can be much different especially in the beginning of the 
outburst, where the spectral modelling indicates that -Rm ^ Rco 
and therefore ojk < tj,. In this case, there should be a region 
around -Rm where a large jump fr om the Keplerian rota tion to the 
co-rotation should occur (LRB95. iLovelace et al.ll2010n . An esti- 
mate of the size of this "transition region" Ar was given by LRB95 

7c^^. (9) 



Ar 



l^Kl^t 



For a thin accretion disc Cs/vk ~ H/r <^ 1, so at r — -Rm the 
transition region Ar ^ H <^ Rm assuming that % ~ i^t and 
7c ~ 1. Therefore, we speculate that its effect on the observed ra- 
diation flux is small due to the small size, but the exact value of 
the correction factor ^ is uncertain. Furthermore, as the accretion 
disc is irradiated by the emission of the hotspot and part of the hard 
emission is absorbed by the disc, the observed colour temperature 
can be altered and therefore affect -Rdisc ■ But because the amplitude 
of Compton reflection SR is poorly constrained, it is not possible to 
accurately take this effect into account. All these unknown factors 
cause a systematic error in -Rdisc, which makes an accurate estima- 
tion of /?m currently infeasible. 

Similarly, the hotspot temperature Thh is also a colour temper- 
ature (at the infinity), but /coi is most like ly the order of unity in 
this case (see lGierlinski & Poutanenll2005l and references therein) 
and the actual hotspot size can be larger d epending on the stellar 
compactness and the hotspot geometry (e.g. Poutanen & Gierlihski| 
^003; Gierlinski & Poutanen 2005; Ibragimov & Poutanen 200^. 
These considerations only change the reported results quantita- 
tively, but the qualitative result of the hotspot size variation shown 
in Fig.[3]is not affected. 

(iii) In our modelling we did not consid er any com plicated spot 
shapes. Based on the simulations of Romanova et al.l 1 12004) a cres- 
cent shaped spot might have been more appropriate to use. How- 
ever, this would have required several new parameters to the model, 
which is not justified based on the quality of the data. This cannot 
be improved because the pulse profile changes so rapidly that co- 
adding more data would not be appropriate. Also by selecting a 
broader range in energies than the 9.8-23.2 keV band to improve 
the statistics cannot be done easily, because the blackbody compo- 
nent (with a different emission pattern) will start contributing from 
below and there are not enough photons above this range. There- 
fore, the only way to address this issue is by making extensive sim- 
ulations and to extract pulse profiles from all the outbursts of SAX 
J1808. 4-3658 and then simultaneously fit them. 



6 SUMMARY 

We have studied spectral and pulse profile variability of SAX 
J1808.4-3658 during its 2008 outburst. We found that the large 
drop in the pulse amplitude and the associated change of the pulse 
profile on September 27 (see fig. 1 in Hartman et al. 2009) was ac- 
companied by a simultaneous spectral transition in the accretion 
disc, which was most evidently seen in the Swift data. Our inter- 
pretation of this transition is that the magnetospheric radius -Rm 
changes because of opening of the field lines that connect the NS 
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magnetic field to the accretion disc. We speculated that the physical 
origin of the field line opening is related to the change of magnetic 
diffusivity rj% or viscosity ut. Through this interpretation we could 
explain why the apparent inner disc radius - which we associate 
with the magnetospheric radius Rm ~ -Rdisc - decreases from 50 
to 30 km in the spectral modelling, why the accretion disc temper- 
ature increases from 0.3 to 0.4 keV and why the disc contribution 
to the bolometric flux increases from 30 to 40 per cent. 

We also saw a decrease in the apparent radius of the hotspot 
in the spectral analysis and this is also consistent with our interpre- 
tation. The field lines that connect the disc at Rm are the ones that 
connect closest to the magnetic pole of the NS. As we assumed that 
the hotspot shape is a ring around the magnetic pole, the opening of 
the field lines should change the location of the inner hotspot edge 
Pin. Therefore, as the outer hotspot edge pout should not change in 
the transition - because for a dipole magnetic field its mostly deter- 
mined by the mass accretion rate that remained constant - we can 
associate the decrease in the hotspot radius to a change in the value 
of Pin. Our pulse profile modelling showed that this is exactly what 
is required to produce the observed transition. Furthermore, our re- 
sults from the pulse profile modelling for the other unknown param- 
eters, such as 9 and i, were also consistent with previous estimates. 
We did not try to place constrains on the neutron star mass and ra- 
dius based on these data, but the identification that varying hotspot 
size can lead to pulse profile changes has profound implications 
for determination of these most important parameters. However, in 
order to place tight constraints, we should not only concentrate on 
modelling these two specific pulse profiles, but instead make use of 
all the outbursts of SAX J1808.4-3658. A detailed re-analysis of 
these data will the main attention of our future work. 

To summarize, for the first time we have found evidence that 
a sudden pulse profile transition in SAX J1808.4— 3658 was most 
likely caused by a change in the way that the NS magnetic field 
is coupled to the accretion disc. This mechanism can be one of the 
causes for pulse profile variability (and the associated timing noise) 
in other AMPs as well. We constrained the disc truncation radius 
at Rt = 21I3 km and estimated the magnetospheric radius to be 
about i?in ~ 30-50 km. This would allow in principle to estimate 
the torques acting on the NS, but the spectroscopic determination 
of i?ni suffers from several systematic uncertainties that currently 
prevent such an attempt. This, however, could be improved if radial 
profiles of uj, B^ and B^ in the magnetospheric region (_Rt < r < 
Rm) were accurately known. 
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